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Abstract 

We consider a deformed flat-band Hubbard model under a periodic boundary 
condition in arbitrary dimensions. We show that the ground state is only all-spin-up 
or -down state. We obtain upper and lower bounds of the one-magnon spin-wave 
energy with an arbitrary momentum. This dispersion relation is the same as that 
in the XXZ model in the certain parameter region. Therefore the spin wave has a 
finite energy gap. These results suggests the our model and the XXZ model. 

Keywords ferromagnetism, flat-band Hubbard model, exact solution, quantum 
effect, anisotropy effect, spin- wave 

1 Introduction 

The origin of the ferromagnetism has been mystery for a long time. We cannot show 
the ferromagnetism without quantum many body effect since interactions for electrons 
are almost spin independent in the microscopic point of view. Moreover, only a non- 
perturbative analysis can predict the ferromagnetism since perturbative approaches such 
as the Fermi liquid theory show only paramagnetism. 

Hubbard model is one of the simplest models to describe itinerant electrons in solids. 
This model is a lattice electron model which considered only two effects. One is an 
electron hopping between lattice points and the other is an on-site repulsive interaction. 
We identify this interaction with the Coulomb interaction. Recently, Mielke and Tasaki 
independently have proposed the models whose ground state has saturated ferromag- 
netism by a rigorous constructive approach. These models are called flat-band Hubbard 
model m El El El- Some remarkable results for ferromagnetic ground states have been 
obtained in this class of models. Nishino, Goda and Kusakabe have extended their re- 
sult to more general models [5^. Tasaki has proved also the stability of the saturated 
ferromagnetism against a perturbation which bends the electron band [HI EI- Tanaka and 
Ueda have shown the stability of the saturated ferromagnetism in a more complicated 
two-dimensional model in Mielke's class jS]. Tasaki has studied the energy of the spin- 
wave excitations in the fiat-band Hubbard model jH]. He has shown that the dispersion 
of the one-magnon excitation is non-singular in the fiat-band Hubbard model, contrary 
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to the Nagaoka ferromagnetism. The flat-band ferromagnetism is beheved to be stable 
against a small perturbation or change the electron number density [10] . 

We argue an anisotropy in a ferromagnetic Hubbard model in this paper. A realistic 
ferro magnet has anisotropics which have no SU(2) symmetry. The Hubbard model has 
SU(2) symmetry and therefore it has no anisotropy. One of the easiest ways to introduce 
an anisotropy is adding a spin anisotropy term such as the Ising term. However, this way 
gives us no information for the origin of the anisotropics as well as the ferromagnetism. It 
is believed that the anisotropy is originated from non-trivial electron dynamics. Moriya 
has shown that the spin-orbit coupling gives a spin-dependence of the electron hopping 
and the effective spin model has the Ising and the Dzyaloshinski-Moriya interactions. 
The hopping anisotropy induces the anisotropy in an effective spin model. We expect 
that the spin-orbit coupling induces the anisotropy even for strong ferromagnetic systems. 
Recently, a deformed fiat-band Hubbard model with an exact domain wall ground state 
was proposed fH^l^j. The deformed model has a hopping anisotropy which depends 
on a spin of an electron. The ground states of the model with open boundary condition 
has same degeneracy as that in original SU(2) symmetric model. The anisotropy induces 
the deformation of SU(2) spin algebra, and then a ground state has a domain wall. It 
is also proven that there exists local gapless excitation above the domain wall ground 
state jlSj. These properties of deformed fiat-band Hubbard model are very similar to 
those in the XXZ model with critical boundary field. Therefore, we expect that the Ising 
anisotropy in quantum spin models comes from hopping anisotropy in a ferromagnetic 
Hubbard model. 

In this paper, we study a deformed fiat-band Hubbard model under a periodic bound- 
ary condition. We can construct all exact ground states. We find that the magnetization 
of ground states are fully polarized, i.e., the ground states are only two states: all-spin-up 
state and -down state. We also obtain upper and lower bounds of one-magnon spin-wave 
excitation energy. We employ the method proposed by Tasaki. He showed that the spin- 
wave excitation in the SU(2) invariant Tasaki model has the ordinary spin-wave dispersion 
relation which has no energy gap above the ground state. We find that the one-magnon 
spin-wave energy has the finite gap above the all-up state in our deformed model. The 
dispersion relation is the same as that in XXZ model in the certain parameter region. 
These facts indicate that our model is related to XXZ model which is non-singular spin 
model. 

This paper is organized as follows. In section [21 we define a deformed fiat-band Hub- 
bard model and show the main results which consist of three theorems. In section El we 
prove the first theorem for ground states. In section |31 we show two lemmas and prove the 
second theorem for the lower bound of spin-wave energy from these lemmas. In section 
and El we prove the lemmas. In section [71 we prove the third theorem for the upper 
bound of the spin-wave excitation. 

2 Definitions and Main Results 

In this section, we define the (i- dimensional deformed fiat-band Hubbard model. We also 
show our main results and discuss their physical meanings. The proofs of results are given 
in later sections. 
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2.1 Lattice 

The lattice A on which defined our deformed Hubbard model is decomposed into two 
sublattices 

A = U A'. (1) 
Ao is a (i-dimensional integer lattice with linear size L defined by 



Ao := |x = G \x, I < J = 1, 2, ■ ■ ■ , d| 



(2) 



where {xM^^ := (xi, 0:2, ■ ■ ■ , Xd)- A' can be further decomposed to Aj (j = 1, 2, ■ ■ ■ , ci), 



i.e. 

d 



A'=Ua,- (3) 

Aj is obtained as a half-integer translation of A^ to j-th direction, 

:= [x + e^^^\x e K] (4) 

where e'--'^ is defined by 

e(^^=(H-,0f=i = (O,---,O,i,O,---,O). (5) 

T 

j-th 

We show the two-dimensional lattice in Fig. [T] as an example. 




Figure 1: The two dimensional lattice (with L = 3). The white circles are sites in Aq and 
the black dots are sites in A'. Electrons at a site can hop to another site if this site is 
connected to the original site with a line or a curve. 



2.2 Electron Operators and the Fock Space 

The creation and annihilation operators for an electron cj, ^ and c^^a obey the standard 
anticommutation relations 

{CxjCT; = 6x,ySa-,T^ {Cx,a: Cy^r} = = {cj, c]^ ,^}, (6) 
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where {A,B} = AB + BA, for x,y & A and electron spin coordinates cr, r =1, i- We 
define tlie no-electron state \E'vac by 

Cx,a'^vac = for Vx e A and cr =t, i • (7) 
We construct a Fock space spanned by a basis 



\.xeA / \xeB 



A,B CA 



We also define a number operator n^. o- by n^^a = cj. ^c^^a whose eigenvalue represents a 
number of electrons at site x with spin a. Note anticommutation relations {cj. ^, cj. ^} = 
z.e. cj.^cj.^ = 0. This relation implies the Pauli principle. We employ the periodic 
boundary condition. This is realized by c^+2Le(i\a = (^x,a for j = 1, 2, ■ ■ ■ , 



2.3 Deformed Flat-Band Hubbard Model 

The Hubbard model is a model which represents a many-electron system on an arbitrary 
lattice. Here, we define a d-dimensional deformed fiat-band Hubbard model. Our model 
is a generalization of the Tasaki model given in the reference |1].A generalized Hubbard 
Hamiltonian consists of two terms 

H := iJhop + -f^int- (9) 
The hopping term ifhop is defined by 

d 

o-=T,i i=i ^-'SAj 
and the interaction term Hi^t is defined by 

Hint = U^n^^^n^^l (11) 

where d^^a is defined on x E Aj {j = 1,2, ■ ■ ■ ,d) 

dx,. := (g^^^'^^/')*c,_,o),, + Ac.,. + (g^^'^)/')* w.),., (12) 

with t,U > 0. q is an complex parameter and p{a) is defined by p{a) = ±1 for a =t, |. 
We write the phase factor of q by 6. The hopping Hamiltonian -f/hop can also be written 
in the following form 

-f^hop = ^ ] ^i,yC^',(TCj/,(T- (13) 

Each term ti^^cj, o-Cj/,. in the hopping Hamiltonian represents the hopping of an electron 

with spin a from site x to site y with a probability proportional to jti'^p. We expect that 
the spin-orbit coupling can induce this spin dependent hopping amplitudes [14 . 

Since the interaction Hamiltonian ifint represents a on-site repulsive interaction, we 
regard this Hamiltonian as a simplification of the Coulomb interaction between two elec- 
trons. 



4 



Note that this system conserves the number of electron. The total electron number 
operator Ne is defined by 

xeK o-=T,i 

Since the Hamiltonian commutes with this operator, we can set the electron number to an 
arbitrary filling. In the present paper, we consider only the case that the electron number 
is equal to |Ao|. We confine ourselves to the Hilbert space Ti spanned by the following 
basis 



n ^-,T n ^-a A, 5 c A with \A\ + \B\ = |Ao| \ . (15) 

Let us discuss the symmetry of the model. An important symmetry is a U(l) symme- 
try. We define spin operators at site x by 

5« := E ct^^^c,., (16) 
where P^'^ (/ = 1,2,3) denote Pauli matrices 

• - C ;) ' - (o \) ■ (-) 

The Hamiltonian commutes with the third component of total spin operator 

[H-, 5'tot] = HS^^^ - sIIIh = 0, (18) 

with 

^S = E^^'- (19) 

Note that this symmetry is enhanced to an SU(2) symmetry in the case of q = 1 i.e. 
Hamiltonian commutes with any component of total spin operator. In this case, this 
model becomes the original fiat-band Hubbard model given by Tasaki |31 Ej. Another 
important symmetry is Z2 symmetry which is generated by a product of a parity and spin 
rotation defined by 

H = H-i = Pexp(^7r^S), (20) 

where P is a parity operator defined by Pc^^aP = c^x,a and Pel ^P = c^-x,cr- n transforms 
Cx,a and cl ^ to c^x,a and cl^.^., where a =| if a =t or a =t if o" =J.. Note the following 
transformation of the total magnetization HS'to/H = —S^^i. An energy eigenstate with 
the total magnetization M is transformed by H into another eigenstate with the total 
magnetization — M, which belongs to the same energy eigenvalue. 

2.4 Ground States 

First, we introduce two states \E'| and \E'| defined by 

= ( n ^It ) ^--^ "^^=(11 ^li ) "^-^ (21) 



where a^^^ with x G Aq is 

a,, = t + ^''t - 5-"""' E <«u,„- (22) 

We can easily verify that the above two states are ground states of the model since the 
Hamiltonian is positive semi-definite and [dy^aiCi\, r] — V ^ ^' ^ ^ ^o- The 
following theorem states that the ground states of the model are only above two states. 

Theorem 2.1 (Ground State of the Model) In the deformed fiat-band Hubbard model 
defined by ^ under a periodic boundary condition with a fixed electron number \Ao\, the 
space of ground states are spanned by two fully polarized states and "ifi- 
The proof is given in the next section. 

This theorem shows that the deformation of electron hopping destroys the large de- 
generacy of ground state in the original SU(2) symmetric Tasaki model and we obtain two 
ground states: all-up state and all-down state. The ground states breaks Z2 symmetry. 
This is similar to that Ising anisotropy effect of the Ising-like XXZ model. 



2.5 Spin- Wave Excitations 



Before we discuss the spin-wave excitation, we remark properties of the spin-wave state 
in quantum spin models. The one-magnon spin- wave state $sw(^) with a wave- number 
k E )C satisfies 

T,^sw{k) = e-^'^-^^swlfc) (23) 

and 

'^S^SW = (^'n.ax - 1)<^'SW, (24) 

where x G and the translation operator is defined by 

TxCy^aT^ ^ = Cx+y,a and Tj-c^y^^T^ ^ = cl._^_y^^. (25) 



K, is the space of wave number vectors 

/C : 



j 27rn 









L-1 L-1 



(26) 



(27) 



2 ' 2 

Then, the one-magnon spin wave state should be in the following Hilbert space Tik 

nk:={^en\ T,.^ = e-'"-^^ and = |(|A„| - 1)^} . 

We define one-magnon spin- wave state as the lowest energy state in Tik- We denote 
the spin- wave energy with wave- number k by i?sw(^)- We can show the following two 
theorems. Note that 6' G M which appears in foUowings is the phase factor of q, i.e. 

I I if) 

q = \q\e . 

Theorem 2.2 (Lower Bound of the Spin-Wave Excitation) There exist positive constants 
ti, Ai, Ui and Ai independent of the system size such that 



Eswik) > 



2U_ 
A^ 



d{\q\ 



-1^ 



J2 cos {2k ■ e(^) + 9) 



A 



(2^ 



for t > ti, A > Ai and U >Ui. 
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Theorem 2.3 (Upper Bound of the Spin-Wave Excitation) There exists positive constant 
A2 such that 



Eswik) < ^ 



(29) 



fort,X,U > 0. 



Remark Theorems 12.21 and 12.31 show that the dispersion relation of the spin-wave is 
asymptotically 



2U 



(30) 



for large A. This dispersion relation is the same as that in the XXZ model with the 
Dzyaloshinski-Moriya interaction given by the Hamiltonian 



+ 



sin 9D ■ (Sr X 5, 



(31) 

with an exchange parameter J = 2f//A^, an Ising anisotropy parameter A = (Iq'I + 
|g|~^)/2 and a Dzyaloshinski-Moriya vector D = (0, 0, 1). Furthermore we can obtain the 
representation (jHT|) as an effective Hamiltonian by a perturbation theory from our model 
with a perturbation parameter 1/A. 

Corollary 2.4 (Existence of the Spin- Wave Gap) There exists positive constant A3 such 
that 



mm Eswik) > 



keJC 



2U_ 
A^ 



di\q\ + \q\-^ - 2) A, 
2 A 



> 



(32) 



This means that there is a finite energy gap between the fully polarized ground state 
and a spin flipped state. We expect that the energy spectra have finite gap above the 
ground states. On the contrary, there exists a gapless excitation in the model under an 
open boundary condition. This drastic difference is also similar to the XXZ model. 



3 Ground States 

In this section, we obtain ground states with the fixed electron number N^. 
prove the Theorem 12.11 on the basis of Tasaki's construction method |4j. 



|Ao|, and 



3.1 Localized Electron Operators 

First, We introduce localized electron operators, which are convenient to construct a 
ground state and to prove the bounds of spin-wave excitation energy. This representation 
was introduced by Tasaki[9j to construct the bases of single electron state. We show the 
construction of the operators in AppendijXl 
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First we introduce a localized electron operator ^ defined by 



(33) 



where is represented by 



-p(a)/4 d 



-o^_„(3).„ + Oi'.j/ H ^ 0. 



^S^+eij),y forxGAo 



A 



A 



x+etj) ,1/ 



(34) 



for a; G A, 



o-^E'vac with X G Ao is a ground state of single electron state. The single-electron ground 
state has |Ao|-fold degeneracy. This macroscopic degeneracy of single-electron ground 
state is one of the origin of flat-band ferromagnetism. {aj. ^\E'vac}a;GAo is a basis of zero- 
energy single-electron states, and {fly^o-^vaclyGA' is a basis of excited single-electron states. 
Furthermore, the ground state energy is and the lowest excitation energy eigenvalue is 
tX^ (See AppendixEl). 

We can also construct another localized electron operator bx^a which satisfies 



{bcc,a,al .^} = Scc,ySa,r aud {&x,<7, fej/,r} = = {al^^,al .^}. 
We represent b^^a in terms of the original electron operator by 

yeA 

dni, (ESI) and (jSni) means 



(35) 



(36) 



(37) 



weA 



Note, ijjy^j cannot be represented closed form but decay exponentially as — becomes 
large, where \\x — y\\i := J2'i=i ^ Vjl- We can obtain the following bounds. 



Lemma 3.1 There exists positive constant Bi such that 

\^ _ I < ^ 



Hi 



< 



oGA 



51 

A2 



(38) 
(39) 



jor all x,y G A. 

Lemma 3.2 There exist positive constants B2, -B3 and B4 

2d 



,(0) 



_A2 + d(|g|l/2+ |g|~l/2)_ 



\\y-^\\i 



for Vx, y G A, 



(40) 
(41) 
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and 



I < 



A3' 



(42) 



The proofs of these lemmas are given in IA.41 

The original electron operators can be written in terms of a|, ^ and bx^a, 



,. = E(^S)*«L and c,,. = 5^^(^6,,.. (43) 

yeA j/GA 



The Hilbert space with |Ao| electrons is also spanned by the basis 



n<r n<ip-- 

,xgA / \xGB / 

because cl „ can be written in terms of at 



A,B CA with \A\ + \B\ = \Ao\ } , (44) 



3.2 Proof of Theorem ED 

Here, we construct ground states of the model. Let \E' be a ground state with |Ao| electrons. 
First, we expand a ground state \[' into the following series 

^ = j2^{A,B) ( n<T ) ( n4.i ) ^vac, (45) 

A,B \xeA / VyG-B / 

where the summation is taken over all A,B G A with \A\ + \B\ = \Ao\. Note that a^^ 
for a; G Ao creates an electron which has the lowest energy of H^op- For the lowest energy 
state of -ffhop, V'(^) B) should vanish, if A or 5 contains a site in A'. Next, we consider 
the interaction Hamiltonian. If we find a state \1/ such that Cx^iCx^i"^ = for Vx G A, then 
the state is a ground state of -ffmt- For x G Aq, we find that ipi^yB) survives only for 
An B = ^, if A, B G Aq. These facts allow us to represent \E' the following form: 

^ = E<^(^)( IK-J'^vac (46) 

o- \xeAo / 

where the summation is taken over all possible spin configurations a = {(Tx)x£Ao- To 
satisfy the condition Cy^iCy^i"^ = for y E Aj (j = 1, 2, ■ ■ ■ ,d), the coefficient holds 

0(^) = g['^'^-=<^'^"'^'^^-'^'^]^V(S-eU),,+e(.)), (47) 

where ax,y is spin configuration obtained by the exchange cr^ and ay in the original con- 
figuration a. The periodic boundary condition allows no configuration which satisfies the 
condition (jlTj) except in the two cases: ax =t for all x G Aq or cr^ =| for all x G Ao- Thus, 
we conclude that all ground states in the periodic boundary condition are only two fully 
polarized states 
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and 



^= n <i I (49) 



This fact proves Theorem 12.11 I 



4 Lower bound of the Spin- Wave Excitation 

In this section, we show that the one-magnon spin-wave excitation has the same dispersion 
relation as that in the XXZ modeL Our proof is based on Tasaki's argument for the 
SU(2) invariant model P|. He has proved that the one-magnon spin- wave excitation in 
the Tasaki model has the same dispersion relation as that in the ferromagnetic Heisenberg 
model. These spin-wave excitations in both models have no energy gap, since they are 
the Goldstone mode above the ground state which spontaneously breaks the SU(2) spin 
rotation symmetry. On the contrary, we will find an energy gap in our anisotropic model. 



4.1 Another Hopping Hamiltonian 

To estimate a lower bound, it is convenient to introduce a new hopping Hamiltonian -ffhop 
which satisfies 

^hop4,a^vac = for X e Ao (50) 

and 

^hop4,,^vac = tA^al^^^vac for X e A'. (51) 

Note that aj. g.^'vac with a; G Aq is the basis of single electron ground states and ,^\l/vac 
with e A' is a basis of single electron excited states. Moreover, the ground state energy 
is and the lowest excitation energy eigenvalue is tX"^, i.e. H^opCil cr^vac = for all a; e Aq 
and (\l/, ifhop^)/||^||^ > tX^ with \l/ = ^ygyv' ^y'^y,cr^vac for any set of complex numbers 
{Cy}yi=\'. Then, we obtain that -ffhop < -f^hop- -^hop is represented by 

Note that H^op is not hermitian and therefore its eigenvectors are not orthogonal systems. 
Nevertheless, all eigenvalues of Hhop are real. We define a new Hamiltonian H := -ffhop + 

We use a representation of interaction Hamiltonian in terms of the localized electron 
operators to evaluate the bounds 

-f^int = Uxi,x2;X3,X4(^x-^^l^CLl2,lbx3,lbx4,1, (53) 

a;i,a:2,X3,X4£A 

where Uxi,x2;x3,x4, is defined by 

Ux,,x2;x3,x4 := u E(v^i:T^v^i:f)*^avi:f (54) 

weA 
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4.2 Basis for the Spin- Wave State 

To define convenient basis of Hk, we define a state for = 0, 1, ■ ■ ■ , d and for a 

set A C A with \A\ = \Ao\ - 1 by 



(55) 



where e^^) = o = (0, 0, ■ ■ ■ , 0) for /i = and e^^) = e^^^ for fi = j {j = 1,2,- ■ ■ ,d). This 
state satisfies both properties and (j21|)- We define another state Q{k) by 

n{k) = ^Y1 ^''""T^^li^oA'fm ^o,AAW, (56) 

which is an approximation of the spin wave state. We will choose a positive constant a 
in the proof. We define the following basis of Tik by 

Bk := Mk)} U i^^^Ak) \ fi = 0,l,--- ,d, Ac A 

with \A\ = \Ao\ - 1 and (/i, A) ^ (0,Ao\{o})}. (57) 

4.2.1 Basic Lemma for Lower Bound 

We define matrix elements /i[\E', $] between \E', $ G by the unique expansion 

M = (58) 

And we define -D[$] by 

D[$]:=3?[M$,<|.]]- 5^ IM-^",^]!- (59) 

We can prove the following lemma. The proof is given in Appendi^Bl 

Lemma 4.1 Let Eo{k) be the lowest energy eigenvalue of H in the Hilhert space Tik- 
Then, we have 

^o(fc) > min (60) 



4.3 Proof of Theorem IZ^ 



First, we show the results of estimates. We obtain these bounds from the direct evaluation 
(See Section El and ini) 

Lemma 4.2 There exist positive constants Fi, F2, F^, F4, F^ and Fq such that 



2U 



A4 



D[n{k)] > 



d{\q\ 



\(1 



^cos(A; ■ 2e(^) + 









F2 


FA_ 


UF^ 


X 


AV 


a A4 


F, 


UF, 




A 


a A2' 





F\ 
A2 



Ua- 



for fi^O and An A' ^ 0. 

We give the proof of this lemma in subsection [HI 



(61) 

(62) 
(63) 
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Proof of Theorem 12.21 from Lemma 14.11 and 14.21 From these inequalities and good 
choice of a, min^gg^, D[^] is given by D[Q{k)]. If we choose a = f//A^, then we find 



F2 + F3/A 
A 



and 



D[^,,a] > tAM 1 



A 



If A > A2 = 4(F2 + F3/A2) and U > Ui = 4F4, then 
If t > ti = 4F6 and A > A3 = [UF^/iAt)]^/^, then 



tA^ 



(64) 



(65) 



(66) 



(67) 



If A > A4 = max {[d{\q\ + \q\-^ + 2)]l/^ [dU{\q\ + \q\-^ + 2)/t]i/6}, then 

m<^i\Q\ + k\-'+^)<^ 



and 



D[n] < 



A4 



A^' 



q\ + |gr^ + 2) < 



2 
tA^ 



We set Ai = max{A2, A3, A4}. Thus we obtain 

2U 



mm Dm = D[Q{k)] > 

'S'&Bk A^ 



(68) 
(69) 

(70) 



for X > Xi, U > Ui and t >ti where Ai = F2 + Fi/Xi. Lemma f4. II and this show that 



Eoik) > '-^ 



^^M±M:!)_^eos(2fc.e(^) + ^)-4l 



(71) 



Since H > H, then Esw{k) > Eo^k). Thus we obtain 

d 



i?sw(/c) > ^ 



^M^^-X:cos(2fc.e(^) + ^)-^ 
i=i 



(72) 



This concludes Theorem 12.21 I 



5 Calculation of Matrix Elements 

In this section, we calculate the matrix elements defined by ()58p to estimate -D[$] for 
$ G Before showing the details of calculation, we list the summary formulae of 
matrix elements. We often abbreviate the fc-dependence of the states for simplicity. 
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5.1 Summary of Matrix Elements 

Here, we show only summary of the calculations. We define hhop[^, ^] and /iint[$, ^] by 

^hop^= J2 ^hop[$,^]$ (73) 

and 

^int^= Yl ^int[$,^]$. (74) 

It is convenient to define and '^i_i,y^v,w with x,v,w G Aq, v ^ w and ?/ G A' by 
and 

^/.,y,i,,«,(fc) := XI e''''""^Li+e(M)a^^+'"i'T^«'+«'i,T4+«'i,T^T oc ^m,(AA{^,«;})uM- (76) 

uiiGAo 

Only and contribute to matrix elements related to /iint[^, '^'] and /lintf'^', 

We summarize the calculations of the matrix elements. 

Vp[^p,B, = ty^WA n A'l + 7^ 0])5^,,x[^ = 5], (77) 

^int =6 ^ ^ ^6 ^f^t,,e(p);j/-x+e(A'),D ~ ^ ^D,e(p) ;t)-x+e('') ; (78) 

and 

where the indicator function % is defined by x^rue] = 1 and x [false] = 0. 

5.2 Treatment of the Hopping Hamiltonian 

It is easy to calculate matrix elements of hopping Hamiltonian, since the modified hopping 
Hamiltonian if hop is "diagonalized" in terms of localized electron operators and 
We can easily obtain the matrix element (f77j) from 

^hop^M,A = t\^{\A n A'l + ^ 0])^^,A. (81) 

5.3 Treatment of the Interaction Hamiltonian 

Here, we calculate the matrix elements. 
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Calculation of /iint['^'; ^/^.x] First, we consider ifint\l/^ .j., 



«'2GA ?J,Ul4gAo 



ik-v 



E E E 

iiii€A' iii26A ii,t04SAo 



(82) 



Here, we decompose W2 into w G Aq and e^^\ We make the change of variables = w + y 
and f = w — X + f for first line. Then we obtain 



E E 

p=0 v,w,y^\o 



ik-{v—x+w)TT t /, 



ik-{w+v—x)TT 

^D+w,«)+e(p);D-x+«;+e(^'),io+i/"^_|_e(p),J,'^"'+J/:T 

ik-(v—x) 



p=0 v,yeAo 
p=0 v,yeAo 



^ik-{y—x)TT ik-{v—x)TT 

^ v,e(p');y-x+e(t'\v ^ ^ v,eM ■y-x+e'-'^') ,y 



^3) 



where we used the relation U. 



x+u,y+u;v+u,w+u ^x,y 



Uxyvw with u G Ao. Thus we find the 



representation (fTHj). I 

For second line of (jH2I), we make the change of variables W2 = w' + e'^p\ wi = w' -\- y, 
Wi = w' + w and v = w' — x + v. Then we obtain 



ik-{v+w'-x)fT y 
/ , / J / J ^ "-^ ?/+«)', -uj'+etP) ;D+t«'-a:+e(^') , to'+i« ^ 

p=0 yeA' ti,to',uieAo 
d 

P=o j/eA' v,weAo w'eAo 

d 

p=0 wgA' v,y&Ao 



?4) 



Thus we find representation (fTUj) . I 
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Calculation of h[$p^y, '^fj,^x,v,w] Next, we consider Hi^t'^^ . 



^i2,i^"'3,i'^l,g+e(M),|'^!oi,T^"'5+1«,T^«'5+f,T^"'4,T'^!u5+I,t'^T 

d 

~ ^ ] ^ ] ^ ^ ^ ^^«)i,'!i;2+e('');«>5+e(f'),ii;5+a::^ 



p=0 UIlSA ti)2,«'5€:Ao 

«)2+e(''),i 

d 



a\ , a|„^^^6^5+u,,-f&u,5+t,,T'^'T + other terms 



p=0 ?i)2,M'5GAo ^ 

,w2+eM -W5;e(p'> ,xb-w5+w,i ) ^'t + other terms 



/ J ^ "^2+e(p),J, 1 '^w,«'5+e(P);e(M),3::"i«2+W5+'U,T 

p=0 W2,W5&ho ^ 

^«;,«)5+e(p);e(M),x^«;2+i"5+»i',T^ "^T + o^her terms 

6 ( ^«),ii;5+e(p);e(f'),a;^p,tU5+''~ 

p=0 uisSAo ^ 

f^«,«;5+e('');e(M),x*P,«;5+«; ) + ot^er terms 



J2 J2 

Uy^y-u,+e(p);e(f-),x^ ^P,y + Oilier terms, (85) 



p=0 yeAo 



where "other terms" do not contribute to the matrix elements /iint[^p,j/5 ^p.x.u.io]- Thus 
we find the representation (jHUI). I 



6 Estimates of the Matrix Elements 

In this section, we show the estimates of matrix elements and to prove Lemma Before 
showing the details of estimates, we list the results of estimates and give the proof of 
Lemma 14.21 

6.1 Summary of Estimates 

Here, we obtain following bounds: 



2U 

nhint[n,n]] >— 



A2 

i=i 



(86) 
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A2 



^7) 



a A* 



< 



a A2 



p=0 viGAo 



9i 

A 



p=0 2/1 G A' iiijtoiGAo 



A 



(89) 
(90) 

(91) 



and 



3f^[/iint[^M,A,^/.,A]]- 
d 

5^5^X[(P,5) 7^ (/i,A)]|/li,,[vl>^,^,Vl>^,5]| > t/ 

p=0 BCA 



(92) 



where x,v,w G Aq with v w, y & A' and /x = 0, 1, ■ ■ ■ , c?. G/ (/ = 1, 2, ■ ■ ■ , 7) are finite 
constants which do not depend on the system size. 

6.2 Proof of Lemma 14.21 

Here, we show the estimates of each We start to evaluate cases that fi = and 

A C Ao- 

Estimates of D[Q] Since i/hop^ = 0, we can represent 

D[Q]=^h,,,,[Q,Q]]- IhuA^^n 

d 



p=0 viGAo 



p=0 j/iGA' t)i,uiiGAo 



(93) 



Note that n = ^o,o/a and a > 0. Then we obtain (jHH), from eqs. (jEHI), (EHl) and (jHH). I 
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Estimates of D[\l/o,a;] -f^hop^o,x also vanishes, tlien we liave 

d 

D[mo^,] =3fJ[/lint[^0,x,^0,x]] - X[(P,^l) 7^ (0,x)]x[(p,t^l) 7^ (0,o)]|/lint[*0,x,^p,.J|- 

p=o gAo 

d 

]|. (94) 

p=0 '1/1 e A' ui.MJigAo 

Tlius we obtain (jSl from (jHI|), (|HH1), (ffl) and (HH). I 

Estimates of Other Contributions We estimate D[^> ^^a\ with 7^ 0, A fl A' 7^ 
or both. Here, we regard and ^ ^,,y,v,w as ^;,,A„\{a;} and ^p,(Ao\{»),«,})u{y} since the 
differences between them are only sign factor which is irrelevant to absolute value of 
off-diagonal matrix elements. Thus we have 

D[^^^a] =^[Kop[^^.A^ ^p,a] + /^int[^p,A, 

d 

J] J]x[(p,i?) 7^ (0,AAW),(/x,A)]|/ii,4v]>^,^,vl>^,B]| - \hU'^,,AM\- (95) 

p=0 BcA 

Since only ^^y^v,w has non- vanishing off-diagonal matrix elements related VL: hmt\^ ^^Ai ^] 
with A n A' 7^ 0, we have 

hA'^,AM<-^ (96) 

from (j89j) . Then we have 



> - ^G, - (97) 
A a A^ 

from (f77j) and (jHS)), where we use 

/ihop[^p,A, ^m,a] = tA2(x[/i ^ 0] + |v4 n A'l) > (98) 

for the cases 7^ 0, A fl A' 7^ or both. Thus, we find the bounds for all contributions 
and we conclude Lemma [4.21 I 

6.3 Estimates of Matrix Elements 

Here, we estimate the matrix elements. To estimate the matrix elements, it is convenient 
to introduce (p^xl by (p'^xl = ip^,l — iPx)t- Note that <^^x}, = ipx}- for \x — y\ > 1. Lemma ITT] 
and 13. 21 give us the bound of Iv^i^ll and Iv^x^^l- Then, we obtain bounds of the matrix 
elements from representations of each matrix element in terms of (^i^L 
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Estimates of 3fJ[/iint[^I/o,x; ^o,x]] First, we estimate the diagonal part of the matrix ele- 
ments. From the representation (fTSj) . we can write 



hmt['^0,x, '^0,xl _ ]_Sr^ (fj _ Jk-ufr 

JT / J \ ^v,o;o,v fc- '-^ 



u 

EE 



x+v,o\v,x 



From the explicit representation of ipj^ja, we obtain 

hmt[^0,x,'^0,x]_,7(o)7(o}.^ ik-xf7(o) 7(o)- 



u 

d 



A2 



E 



r,(°) 



r.io) 



+ 



A2 



0~l/2 _ -ifc.(x+2eO)) T(o) 7(o) 



1/2 



-x-e(j),T^x+e(j),i 



+ 



ifc.(:j_2e(i)) 7(0) 7(0) 
e '/^_3,.+e(j),T'''^x-e(j),i 



(99) 



(100) 



where we use the relation ^/'i^j'"^ = ip^x-u,a fo^^ ^ ^ ^o- 
For X = o, we obtain 



A4 



A2 



^-1/4 



E' 

i=i 

(g-l/2_gl/2g-.fc-2eO)) + 
^^'''^Z (o) , (o) N , 1 (o) (o) 



A 



X 



^-l/2gifc-2eW) 



) ■ 



(101) 



Lemma EiH means Iv'^^gCj) ^| < -B3A then we find the bound (jHBj) . I 

For X ^ o, we find (|HZD from < 1 - 0(A-2) and = 0{\-^) for x ^ o. I 
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Estimates of /iint[^o,y5 ^0,0] Next, we estimate tlie off-diagonal part of the matrix ele- 
ments. To obtain (jHHjl . we calculate the following matrix element 



U 

A2 



i=i L 



Since iV'i"^] = 0{X~^) for w G A', we find the bound 



(102) 



Estimates of h[i/ fj,^x,v,w, ^0,0] Using the following bound 



U_ 
A2 



< 0{UX~^), 



we obtain (|S^ . 



Estimates of /z[\E'o,x; with x ^ o Using the following bound 



(103) 



J2 Yl l^int[*0,x,^p,y]|x[(p,2/) 7^ (0,X)] 

p=o yeAo 

d 



V T^ui,! i^ji), J, / \^v,x—y'-'w,x—y '-'i;,x+i/'-'tj 



''EE 

x[ip,y) y^iO,x)] +0{UX- 



I^SCl - e-'<^^y^iT^^if^ I x[(P, y) ^ (0, X)] + oiux- 

p=0 yeAo 

U E - e-^'^-^^+^^V^ir^V^S X[(p,y) 7^ (0,x)] + 0(t/A-^) 



< U E (l^^i^a 

I/6A0 

< 0{UX-^). 



+ 0(f/A 



+ 



",(2^)J,(°) 



(104) 
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for X 7^ o, and using another bound 



$^$^|/^int[^0,o,Vl/p,,]|x[(p,l/)7^(0,o)] 

p=o j/eAo 



< 0{UX-^), 



o{ux-^) 



for X = o, we obtain (|9fl|l . 

Estimates of /i[\l/o,x, ^p,y,j;,ti;] Using the following bound 

d 

M=o ?/eA' v,weAo 



j=l yeAj D,uieAo 



j=l yeAj v,weAo 



we obtain (jHH). 



(105) 



(106) 



Estimates of Other Matrix Elements To obtain (|^. we calculate the matrix ele 
ments of Hi^f 



Uli,'(«2i'^3)''i'4SA 



E 

weAo 



n 



EY^ Jk-wfr t , t „t 



uigAo uii,-u;2,W3,ui4gA 



,T^«'4,T n ^i+'^.T 
\veA J 



E e"""U^un.2;w+eM,n.4X[w4 eA + W,Wi^{A + W)\{W4}] X 
ioGAo iiii,ui2,«)4SA 



sgn[A + w,Wi,Wi]al^^^[ JJ 



.v£{{A+w)\{w4})U{wi} 



20 



Here, we make the change of the variables W2 y + e^p\ w ^ w + y, Wi ^ Wi + y and 
W4 W4 + y. Then we obtain 

d 

x[wi + yeA + w + y,wi + y^{A + w + y)\{w4 + y}] x 
sgn[y4 + w + y,wi + y,W4 + y]x 



vei{A+w)\{w4})U{wi}+y 



ik-WTT y 



EE E 

p=0 liigAo «)i,ui4gA 

x[w4 & A + w,Wi ^ {A + w)\{w 4}]sgn[A + w, wi, x 



e ^a^( 



n 



g(p) + !/,j 

yeAo \ fe((A+tu)\{«)4})u{«)i}+y 



~ ^ ] ^ ] ^ ] ^ ^i«i,e(p);tu+e(M),«i4 X 

p=o tueAo «)i,ui4eA 
X[w4 G A + w, wi ^ (A + w)\{w4}]sgn[yl + w, Wi, W4] X 

'^p,iiA+w)\{wi})U{w^}- (107) 

Thus, we find 

(108) 

Here, we neglect the sign of the matrix element, since it is irrelevant to obtain the bounds. 
If we define A' by A' = {{A+w)\{w4})U {wi} , then we obtain its solution A = {A'\{wi})U 
{104} — w. Thus, we obtain 

hint[^p,A, p„(A\{wi})U{w4}-w] = e'''''"U^^^e(p);w+eM,W4X['^^ G A, W4 ^ A\{wi}]. (109) 

If wi = W4 and w = o, then A = {A\{wi}) U {^^4} — w. Since other wi, W4 and w makes 
(v4\{wi}) U {W4} - w = A,we find 

3f?[/iint[^p,A,^p,A]] - Yl l^int[^p,A,$]| 
*eB\{1'p,A} 

|x[(w,ti'i) 7^ (e(''\w4)]. 

.uiisA J «)iGAtoeAi«4e(A\yl)u{tui} 

(110) 

To estimate the first term in (jliop . we use 



3? 



^ ] ^u;i,e('');e(''),-u)i 

,wieA 



0(t/A- 



:iii^ 
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We estimate the second term in ()110|) . From the exphcit form of ip^}-, we obtain 

wi gA weA w4e{A\A)U{wi } 

<UY1 E I^SCi^lx[^4,^i ^ e^'')] + 0{UX-'). (112) 



wieA w4e{A\A)U{wi} 



Since l^'lTV^I — ^(-^ ^) ^^r W4 7^ Wi, then pi2p is bounded by 



< ^ E i^Sli^STi^t^i ^ e^'^] + (113) 



Since iV'i^^Vl — ^("^ ^) ^^"^ bound 

wi g A iogA «)4g (A\A)U{mji } 

<0{UX-^). (114) 
From (ITTT|1 and (fTTll . we find the bound 



7 Proof of Theorem [2T3 



In this section, we show the proof of the upper bound of spin-wave energy. The proof is 
given by the standard variational approach. If we can show 



const 



:ii5) 



this concludes the Theorem 12.31 Here we evaluate this quantity. 
Since i^hop^o.o = 0, we obtain 

(^0,0, ^^"^0,0) 



112 /iint[^0,o,^0,o] + 



0,0 

d 



EE E x[v^w]K^,[^^„,<I^,/^^^^^^ (116) 

Note (\E'o,o, ^j,xO = for j = 1, 2, ■ ■ ■ , d because of the representation of ay^a with y G A' 
in terms of its dual 

y'GA' 

where Cy'^y^^ is a complex coefficient. We also find (^'0,0, ^At,y,t;,u)) = with the same 
argument. Thus, we obtain 

— = /iinti^O.o, *0,oJ + 2^ /imt[^0,a:,^0,oj— ||-r (117) 

x&ko\{o} 
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We already have the bounds of the matrix elements. We concentrate to estimate the inner 
product between \E'o,o and ^E^o.x- We represent (\E'o,o, ^o,a;) 



-ik-{v-w) jij 



}Ki,<i}ll*Tl 



:ii8) 



Using the relations 



«L} = and {bl^, by,.} = J2 ^S(^Sl)* (119) 



and the explicit representation of V'i^CT, we obtain 



A2 



1 



(120) 



To obtain this, we have also used the translational invariance 

ofi^^S-- = with 

1(7 e Aq. Thus we have 



J2 /i[^0,x,^0,o](^0,o,^0,x) 

xeAo\{o} 



^^'tI^3 E l(^o,o,vl>o,.)|<^'^^i^i||^Tll', (121) 



xeAo\{o} 

where we use fj88p and Lemma f3. 11 The norm of \E'o,o can be estimated as 

-A2 + d(|g|l/2+|g|-l/2) 



II ^ 0,0 II 
1 

T2 5Z 



uigA 



:i22) 



Then we have 

(^0,0, i^^o.o) 



0,0 1 



<|/l[vE'o,o,^'o,o]| + ^^ 



<- 



2U 



d{\q\ + \q\ 



-^cos(/t-2e(^) + ^) + 



K1/K2 



(123) 



This completes the proof. I 
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A Localized Basis 



A.l Single Electron Problem 



First, we consider a single-electron problem with periodic boundary condition to construct 
the localized basis. A single electron state with spin a: $0- can be written in a form 



(A.l) 



Since ifint^o- = is obvious, the Schrodinger equation = e^^'^^a- is reduced to 



From the representation of -ffhop- 



xeAo cr=t,J, 



l/2^ 



a 



j = l xGAj a=1,l '- 



(A.3) 



we obtain an expression of the equations 



e^x,a = < 



td{\q\'^' + \q\-'^')^IJx,.+ 

d r 



(A.4) 



if a; G Ao 
if X G A' 



Since 



t^l for w G Ao, we can use the Bloch's theorem which ensures that ipx,a 
can be written in a form 

^x,a = e'^'^Vx^^ik) (A.5) 



where Vx,a{k) satisfies 

for w E Z, and /c G /C is a wave number vector (see eq. 



MX 



if X G Ao 
j if X G Aj 



(A.6) 

We define a mapping u by 

(A.7) 
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We rewrite eq. ()A.5|) to a form 



Ak-x 



(A.8) 



since v^^aik) is translationally invariant. By substituting the representation ()A.8|) into the 
Schrodinger equation ()A.4|) . we find the equation 



(A.9) 



where Va{k) is a + l)-dimensional complex vector defined by Va{k) = (i^/^,(t(/c))^=o 
M^'^^k) = (mS(A;));^^^^o is a 1) X {d + 1) matrix. mS(A;) is defined by 

ik-iy-x) ^^_^Q^ 



for u{x) = 11 and u{y) = p. Here we define Ap = Aq for p = and Ap = Aj for p = j 
for convenience. We use indices p and p as an integer of 0, 1, 2, ■ ■ ■ ,d, and j and I as an 
integer of 1, 2, ■ ■ ■ ,d. 



Mi:}{k) 



rf(|g|i/2 + |g|-V2) + 2 5^cos 



e 



p{a) + 2A; ■ e^^'^ 



A 




/i = p = 

/x = 0,p = j. (A.ll) 

^J' = p = j 

otherwise 



We denote the energy eigenvalues as e^T'' (A;), where we assigned index p such that e^^^ (k) < 



^u+i(^)- We can easily calculate the eigenvalues of the matrix M^"'\k): 



el^\k) 





A2 



A2 + d(|g|V2 + |g|-i/2) + 2^cos 



for yU = 

for /i = 1, 2, ■ ■ ■ , d — 1 







p((t) + kj 



ioT p = d 



We define vf\k) = {vfl{k))%^ by 
1 



v^'lik) 



if /i = 



1 (g-pW/4g-ife.eO) ^ ^pM/4gifc.eO)) ^^^-(^1^2, 



A 



(A.12) 



(A. 13) 



This vector is an eigenvector which belongs to the zero-energy eigenvalue of M^'^\k). We 



also define Va \k) = {vli,a{k))'f^^Q for each j by 



if p = j- 
otherwise 



(A.14) 
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This is orthogonal to v^\k). The vectors {^^ct^''(/^)}^=o fixed k and o" is a basis of C^^"^, 
since v^\k) with // = 0, 1, ■ ■ ■ , d are hnearly independent of each other. 

We also introduce the dual of the basis {v^a\k)}'fj,^Q- We define the Gramm matrix 
GM(A:) = (Gafc))^,^^oby Ha 



G'^:Uk) = {vi^\k)A'\k))- 



(A. 15) 



Since the vectors vi'^\k) with /i = 0, 1, ■ ■ ■ ,d are linearly independent of each other, the 
corresponding Gramm matrix is invertible. We define the dual vectors v^\k) by 



p=0 



for each /x. They satisfy 



(A.16) 



{vi'\k)A'\k)) = 5,,, and = 5,,,. 

r]=0 

A. 2 Construction of the Localized Basis 



We define ijjy^l and ipy^l by 



{x) 



and 



They satisfy 



We can easily obtain explicit representation of ipy^l: 

-p(a)/4 'i 



^ / ^ "x-e^J^y ~^ ^x,y 



-p(ct)/4> 



^ .=1 



(i) „ if z G 



-5 



if X e Kj 



(A.17) 



(A.is: 



(A.19) 



(A.21) 



A.3 Calculation of 

To estimate 'ip'^Jr-, we obtain the explicit representation oiv^\k). We write G*-'^^(A;) in the 
form 

(\v^a\k)\'' ■•■ 0\ 

Q^''\k) = ° 1 

V / 



1 + —w„{k)®w„{ky 



(A.22) 
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where 1 is the d x d identity matrix, and w„{k) and Wfj{k)* are d-dimensional vectors 
defined by 

W,{k) = ('(g-P(-)/4)*e«fc.e(i) ^ (^pM/4)*g-ifc.eO)y ^^ 23) 



and 



(A.24) 



(S> denotes the standard tensor product. There is a general formula to obtain the inverse 
matrix 



-1 



1 



Then, we obtain the explicit form of G'^'^''(A;) 



w,{k)0w^{ky. (A.25) 



G^''\k) ' 



f\vi^\k)\-^ 




1 



k.(A;)|2 + A2 



V 

Thus, we find an explicit representation of vi^\k): 
1 il 



0\ 

Wcr{k) (g) w^{k)* 

J 



1 



A 



for yU = 
for ^ = j 



(A.26) 



(A.27) 



where 



= ^i^i(^) + - — ^^^—^ X 



S, + S!f\k) 

+ ^pM/4g*fc-e(0 /i = /(= 1, 2, ■ ■ ■ , ' 



A2 + 5i + 5j'^)(A;) 
for /i = 



(A.28) 



5, = ci(|g|i/2 + |g|-i/2)^ = 2^cos('2A;-e(^) + -p(a)y (A.29) 

7 = 1 ^ ^ 



A. 4 Proof of Lemma 13.11 

Here we prove Lemma 13.11 We give the estimates of 'i/j^Jt only for x, ?/ G Aq, since the 
estimates for other x, y G A are almost the same as that for x, y G A^. 



From the definition of ip^Jr, we can represent 



Jk-{x-y) _ 



kelC 
1 1 



keK. 



1 _|_ 5l _|_ '^2 ^(^) 



^^ik-{x- 



■y). 



1 + 



(A.30) 
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We can expand the last factor, since \S^''\k)\ <2d<X^ + Si. Then, we obtain 



11 



A2 n=0 keK. 
oo 



Here, 



^ik-{x-y) 



.i=i 



= E 

d 

E 

Thus, we find a bound 

E wsi <r^E (^)" E ^E^-''-"sr'(*)" 

°° / 1 \ " 

EE E 



jl,j2,--- ,in = l Sl,S2,--- ,Sn = ±l 



(A.32) 



^ CXD 

1 4- Si 



2d 



1 J. Si _ 2d- 
A2 A2 



(A.33) 



We can represent ip^}- — ip^x}a- in terms of ipw^a where w G Aq with |w — x| < 1, 



1 I Si _i_ \k) 
^ + 1? A^ 



-2seO),(T' 



(A.34) 



i=i s=±i 



Thus we find the following bound 



< 



const 



(A.35) 



For other x, |/ G A, we can write ^/'y^i in terms of tjjv^a with G Aq like the representation 
^KMf and can obtain the bounds. Then we can conclude the lemma ITTl I 
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A. 5 Proof of Lemma 13.21 

Here, we give proof of Lemma [3 .21 We only show the proof for the case of x,y & Aq, since 
the estimates for other x,y E A are almost the same as that for x,y & Aq. 
In the eq. flA.32|) . if \\y — x\\i > n, then S o,) = 0. Thus we find 

2 ^ / 1 \ ^ ^ 

n=\\x-y\\i ^ ^ ji,j2,--- ,jn=l Si,S2,--- ,Sn=±l 

^ 1 ,22, / 2d y 



^ n=\\x—y\\i 



2d 1 , 

(A.36) 



For other x, y G A, we can find similar bounds. Then, we obtain the first inequality in 
Lemma 13.21 

We also obtain ()41|) and ()42|) from the representations 

fce/c I + j2 ~\ X2 — s=±i 

and 



7(o) /(o) 'S'l ~(o) 1 v-^ isfEM 7(o) /. oo^ 

«=1 s=±l 

They conclude the Lemma f3. 21 I 



B Proof of Lemma 14.11 

Let $0 £ '^fc be the eigenstate which belongs to the lowest energy eigenvalue Eq. We can 
expand $o by the basis Bk- 

$0 = ^ a^^)^. (B.l) 

The eigenvalue equation is 

J2 C{^)h[<l>,^<^ = J2 ^oC{<^)<^. (B.2) 

Thus, we have an equation for the coefficient 

EqC{<I>) = J2 h[^,'^]CW (B.3) 

for each ^ E Bk- Then, we obtain 

*GBfe\{*} 
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Since Eq G M, the real part of the equation is 



cm 



3?[/i[$,<l>]]+ J2 cos(^$,^) |/i[<l>,^ 

C(vl/) 



C($) 



*eBfc\{*} 



cm 



(B.5) 



where ^$^4, G M gives the phase factor of \l/]C(\l/)/C($). If we choose $ such that 
|C(<I>)| = max^eB, |C(^)|, i.e. |C(^)/C($)| < 1 for V* G ^B^, then we find 



> $]] - ^ |/i[<l>, ^] I = 

*eBfc\{*} 

Since -D[$] > min,i,g0^, -D[\I/], we have proved lemma HHJ I 



(B.6) 
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